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ABSTRACT 

We report the observations of the magneto-hydrodynamic waves propagating 
along magnetic flux tubes in the solar photosphere. We identified 20 isolated 
strong peaks (8 peaks for pores and 12 peaks for inter-granular magnetic struc- 
ture) in the power spectra of the l.o.s. (line-of-sight) magnetic flux, the l.o.s. 
velocity, and the intensity for 14 different magnetic concentrations. The obser- 
vation is performed with the spectro-polarimeter of the Solar Optical Telescope 
aboard the Hinode satellite. The oscillation periods are located in 3—6 min 
for the pores and in 4—9 min for the inter-granular magnetic elements. These 
peaks correspond to the magnetic, the velocity, and the intensity fluctuation in 
time domain with r.m.s. (root mean square) amplitudes of 4 — 17 G (0.3 — 1.2 
%), 0.03 — 0.12 km s -1 , and 0.1 — 1%, respectively. Phase differences between 
the l.o.s. magnetic flux (0_b), the l.o.s. velocity ((f> v ), the intensities of the line 
core (07,core)) an d the continuum intensity (0j )CO nt) have striking concentrations 
at around —90° for (fi B — <f> v and <f> v — (pi tCOTC , around 180° for (pi jCorc — 0b, and 
around 10° for 0/ iCO rc — <t>i, cont- Here, for example, 4>b — 4>v ~ —90° means that 
the velocity leads the magnetic field by a quarter of cycle. The observed phase 
relation between the magnetic and the photometric intensity fluctuation would 
not be consistent with that caused by the opacity effect, if the magnetic field 
strength decreases with height along the oblique line of sight. We suggest that 
the observed fluctuations are due to longitudinal (sausage mode) and/or trans- 
verse (kink mode) MHD waves. The observed phase relation between the fluctu- 
ations in the magnetic flux and the velocity is consistent with the superposition 
of the ascending wave and the descending wave reflected at chromosphere/corona 
boundary (standing wave). Even with such reflected waves, the residual upward 
Poynting flux is estimated to be 2.7 x 10 6 erg cm~ 2 s -1 for a case of the kink 
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wave. Seismology of the magnetic flux tubes is possible to obtain various physical 
parameters from the observed period and amplitude of the oscillations. 

Subject headings: Sun: MHD wave — Sun: photosphere — Sun: coronal heating - 
Sun: solar wind 
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1. Introduction 

Alfven waves or more generally transverse magneto-hydrodynamic waves would play 
a key role in coronal heating and solar wind acceleration (e.g. Suzuki & Inutsuka, 2006). 
Numerous studies about generation, propagation and dissipation of the Alfven waves have 
been carried out observationally and theoretically (e.g. Ryutova & Priest, 1993). Alfven 
waves would be generated in the high-/? region of the solar atmosphere. Its precise power 
spectra is, however, not observationally known. Ascending Alfven waves with wavelength 
longer than the Alfvenic scale height may be reflected back at the chromospheric-coronal 
boundary (Moore et al.1991; An et al. 1989; Hollweg 1978; Suzuki & Inutsuka, 2005). It 
is poorly known how much Alfven-wave flux generated in the photosphere is propagated 
all the way to the corona through the fanning-out flux tubes. High-quality observations 
to obtain spectra of magnetic fluctuation is of crucial importance to understand coronal 
heating and acceleration of fast solar wind. 

Ulrich (1996) made the first critical observations, and reported the detection of the 
magneto-hydrodynamic oscillations with properties of the Alfven waves. He suggested 
that the observed phase relation between the magnetic field and the velocity perturbation 
is consistent with the outgoing Alfven waves. The observing aperture of 20" x 20" is, 
however, very large compared with the spatial scale of the flux tubes along which the 
Alfven waves propagate. Such a large aperture may make it difficult to identify the weak 
transverse waves with different frequency and phase, which might become evident in higher 
resolution observations. Velocity and magnetic field oscillations in the sunspot umbra were 
detected by Bellot Rubio et al. (2000), Lites et al. (1998), Norton et al. (1999), Riiedi 
et al. (1998), Riiedi & Solanki (1999), Balthasar (1999), and Settele et al. (2002). Riiedi 
et al. (1998) and Bellot Rubio et al. (2000) obtained the phase difference of -90° and 
90° between the fluctuations of the line-of-sight velocity and the magnetic field strength 
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4>v — 4>b, respectively. They suggested that the magnetic field fluctuation is caused by the 
opacity fluctuations that move upward and downward the region where the spectral line 
profiles are sensitive to magnetic fields. Norton et al. (2001) obtained the center-to-limb 
dependence of the phase angle between the magnetic and the velocity fluctuations with the 
Michelson Doppler Imager aboard the SOHO satellite. They reported that the phase angle 
is near —90° at the disk center and near 0° at the limb, and made an important comment 
that the Alfven waves be more easily observed at the limb. They suggested that the phase 
relation reported in the paper is not due to the opacity effect. Khomenko et al. (2003) 
compared the analytical solution of the MHD equations including gravity, inclination of 
magnetic field, and effects of nonadiabaticity with the observations reported by Bellot 
Rubio (2000), and concluded that the detected time variation in field strength could be 
partly due to magnetoacoustic waves. Riiedi & Cally (2003) suggested that most of the 
observed magnetic field oscillations is due to the opacity effect caused by temperature and 
density fluctuations associated with magnetoacoustic waves. 

Recently apparent transverse oscillations, which are clear evidence of the Alfven waves, 
are detected in prominence (Okamoto et al. 2007), in spicules (de Pontieu et al. 2007, He 
et al. 2009), and in Ca jet (Nishizuka et al., 2008) with the Solar Optical Telescope (SOT; 
Tsuneta et al. 2008a; Suematsu et al. 2008; Ichimoto et al. 2008; Shimizu et al. 2008) 
aboard the Hinode satellite (Kosugi et al., 2007). These Alfven waves have enough Poynting 
flux to potentially heat the corona. We, however, cannot rule out the possibility that these 
waves are the standing Alfven waves. Transverse oscillations of coronal loops are detected 
by Taroyan et al. (2008), Mariska et al. (2008), and Van Doorsselaere et al. (2008) using 
the EUV Imaging Spectrometer (Calhane et al., 2007) aboard the Hinode satellite as well. 
Ubiquitous upward Alfven waves in the corona are detected by Tomczyk et al. (2007) using 
the Coronal Multi-Channel Polarimeter without magnetic field information. We stress 
that the observations of the magnetic field fluctuation with the simultaneous velocity and 
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photometric measurement allows us to identify propagating hydromagnetic waves. 

The literatures so far introduced are mainly concerned with the pure Alfven waves. 
The magnetic fields in the solar atmosphere have a form of isolated magnetic flux tubes 
embedded in a nearly field-free fluid. Such flux tubes carry the incompressible torsional 
Alfven waves, and the linearly polarized Alfven waves can exist only in the uniform 
media. The flux tubes also carry the kink waves (transverse waves) and the sausage 
waves (longitudinal waves) (e.g. Stix, 2002) instead of the linearly-polarized Alfven waves. 
Magnetic tension force of the flux tube is the restoring force in the kink mode (e.g. Spruit, 
1981), and is essentially incompressible. The sausage mode with the azimuthal wave number 
m = 0, as first defined by Defouw (1976) and discussed e.g. in Roberts and Webb (1978) 
and Ryutova (1981), is related to a slow magnetoacoustic mode. In the sausage- shaped 
perturbed boundary of the flux tube, where the flux-tube area increases, the magnetic field 
decreases, whereas the plasma pressure increases; vice versa. A fast magneto-acoustic mode 
propagates across the flux tube, and is not localized radially in the vicinity of the flux tube; 
we do not regard this as a mode of flux tube oscillations. In this paper, we report a clear 
detection of magnetic, velocity and photometric oscillations of the magnetic flux tubes with 
the Spectro-Polarimeter (SP) of SOT. The data is extensively analyzed in terms of both 
the linearly-polarized kink waves and the slow sausage waves, while we will not discuss the 
torsional Alfven waves due to our constraint in the analysis as we explain later. 

SOT/SP is ideally suited to detect the magneto- hydrodynamic waves propagated along 
the flux tubes due to its high spatial and time resolution and its high polarimetric and 
photometric precision (e.g. Ploner & Solanki, 1997). SOT/SP obtains two spectra of iron 
lines (Fe I) with wavelengths of 630.15 nm and 630.25 nm, which are suitable for observing 
lower photosphere (del Toro Iniesta, 2003). Earlier studies about magnetic fluctuations 
were done in sunspot umbra, since small-scale flux tube (~ 1") fluctuations might be 
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difficult to detect. The high spatial resolution of Hinode (~ 0".16) allows us to detect the 
fluctuations in such small-scale flux tubes. Furthermore, stable observations from space 
allow us to detect clear intensity fluctuations for the first time, and to obtain the phase 
relations among the fluctuations in the magnetic flux, the velocity, and the intensity. This 
allows us to examine the opacity effect more in detail. 

For the detection of weak magnetic fluctuations, we prefer to use the Stokes V signal 
instead of the Stokes Q or U signal because of its much higher sensitivity to magnetic flux. 
We, thus, intentionally choose magnetic concentrations located away from the disk center 
to observe possible fluctuation of the transverse magnetic field in the Stokes V signal and 
the associated velocity signal. We estimate the magnetic field fluctuation associated with 
the Alfven waves to be about 5B = 10 G by substituting typical values for the photosphere 
(magnetic field strength of a flux tube B = 2000 G, Alfven speed va = 20 km s -1 , and 
velocity fluctuation Sv = 0.1 km s _1 ) to the relation about the Alfven wave 5v/va = SB/B 
(e.g. Priest 1981). A detection limit of the longitudinal and transverse magnetic fields 
observed by the SOT is known to be 1-5 G and 30-50 G, respectively (Tsuneta et al. 
2008a). This exercise demonstrates that the SOT/SP can detect the transverse MHD waves 
in the Stokes V signal with high signal-to-noise ratio, if such MHD waves are present in the 
photosphere. 



2. Observations and Data Analysis 

2.1. Hinode Observation 

The data used in this paper was taken on 2007 February 3 — 5. The region that 
we observed was NOAA 10940, which moved from west 25.2 to 49.0 degrees in longitude 
during the course of the observation. The region consists of pores and magnetic flux 
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concentrations located in the inter-granular lanes (Ishikawa et al., 2007), which we hereafter 
call Inter-granular Magnetic Structure (IMS). The integration time is 4.8 s, and the field of 
view is 1".92 (EW) by 81". 92 (NS). The pixel size is 0".16. Periodic scanning was done by 
the SOT/SP for 1 hour or 3 hours depending on the flux tubes with cadence of 67 s. This 
time resolution allows us to detect magnetohydrodynamic waves with a period longer than 
134 s according to the Nyquist criteria. 

We analyzed 14 magnetic flux concentrations as tabulated in Table [TJ All these 
magnetic flux concentrations are of positive polarity (magnetic field vector toward the 
observer). The region #05 is shown in Fig. [I] as an example of the data. The region #05 
contains a pore in a plage region. 

2.2. Time-Profile Data Analysis 

We use the Stokes I and V profiles of the Fe I 630.25 nm line to derive the line-of-sight 
velocity, the line-of-sight magnetic flux, and the intensity. The line-of-sight velocity 
fluctuation (5v\ os ) is derived by measuring the Stokes V zero cross position A c . The Stokes 
V profiles reflect the motion of the magnetic atmosphere better than the Stokes I profiles, 
which also contain the information of the non-magnetic atmosphere. The line-of-sight 
magnetic flux fluctuation (5$i os ) is derived with the help of weak field approximation (Landi 
degl'Innocenti & Landolfi 2004) rather than relying on the standard Milne-Eddington 
inversion (e.g. del Toro Iniesta 2003). The Milne-Eddington least-squares fit is performed to 
the observed Stokes profiles of the Fe I 630.15 nm and Fe I 630.25 nm with 12 parameters, 
which may be subject to noise that impedes the detection of fluctuation with amplitude as 
small as SB/B ~ 0.4 %. In the weak field approximation, the line-of-sight magnetic flux is 
proportional to the degree of the circular polarization CP defined by: 



where V(A) is the Stokes profile observed with the SOT/SP, A c is the measured zero cross 
position of the observed Stokes-V profiles as described above, d\ is 43.2 pm, and / con t is the 
continuum intensity. The observed Stokes I and V profiles for the region #5 (Table [1]) are 
shown in Fig. [2] as an example. Since the integration is done with respect to A c , and the 
integration range is wide enough to encompass the entire profiles, the integral should not 
have any cross-talk with the velocity. Intensity fluctuations in the line core (SI core ) and in 
the continuum (5/ cont ) are derived from the line core intensity J core and continuum intensity 
7 cont defined by, 



where /(A) is the Stokes I profile observed with the SOT/SP, d 2 and d 3 are 54.0 pm and 
10.8 pm, respectively, and the factor of 4 is to adjust the difference in the integration range 
between CP and /. 

The intrinsic magnetic field strength (B Q ) and the filling factor / are derived from 
the Milne-Eddington inversion to accurately determine the Alfven speed. The intrinsic 
magnetic field strength B is used only for this purpose. The filling factor is defined as the 
fraction of area occupied with the magnetic field in a pixel (Orozco Suarez et al., 2007). The 
12 free parameters are intrinsic field strength (-Bo), inclination and azimuth for magnetic 
field vector, line strength, Doppler width, damping factor, Doppler velocity, source function, 
source gradient, macro turbulence, filling factor (stray-light factor), and the Doppler shift 
of the stray-light profile. 

We should track the region of interest (ROI), for which the wave analysis is performed, 
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in a Lagrangian way for an extended period of time. In the case of pores, the overall 
magnetic structure is maintained over 1 hour as shown is Fig. [TJ In this case, we set the 
ROI to cover a portion of a pore. The size of the ROIs for pores is typically 2"x 2"to 2"x 
5". The physical parameters are averaged inside the ROI. On the other hand, IMSs are 
generally not maintained for 1 hour: magnetic elements may combine, split, or decay within 
a time period of several tens of minutes. Thus, we set the ROI in this case large enough to 
encompass the entire magnetic flux concentration in the spatial and temporal domain, and 
average the physical parameters of the pixels with CP larger than 0.01 inside the ROI. The 
size of the ROIs for IMSs is typically l"x l"to 2"x 4". Examples of SP images for a pore 
(ID #05) and an IMS (ID #10) are shown in Fig. [3j The average line-of-sight magnetic 
flux $ los is given by 

EN 

where N is the number of pixels inside the ROI for which CP is larger than 0.01. A is the 
conversion coefficient for converting CP to magnetic flux. A is estimated to be 4.16 x 10 -5 
G _1 in the Appendix. Positive values for the line-of-sight magnetic flux and velocity 
indicate that the are directed towards the observer. 



3. Power Spectra and Phase Relation 

The top panels of Fig. [5] show the time profiles of the line-of-sight magnetic flux, the 
line-of-sight velocity, and the line core intensity for the region #04 (Tabled]). We applied 
the Fourier Transform to all time profiles. The result for the region #04 is shown in the 
bottom panels of Fig. HI The power spectra generally show one or two isolated sharp peaks 
in the shorter periods, while broader peaks are found in the longer periods, corresponding 
to a gradual rise and fall in the time profiles. Some of the peaks have the same period in 
the magnetic and velocity field, and the photometric intensity. We found 20 such common 




Fig. 1.— Top: The SOT filtergraph (FG) image taken in the Fel 630.25 nm line at 19:45(UT) 
on 2007 February 3. The field of view is 217". 1 (EW) x 108".5 (NS). The pixel size is 0".108. 
Exposure time is 90 ms. The black rectangular box indicates the region ^05 (Table [1]). 
Bottom: Zoomed SP images (Stokes I and V) for the region ^05 taken at 19:15-20:18 (UT) 
on 2007 February 3. Periodic scanning was done by SP for about 1 hour with a cadence of 
67 s. The integration time is 4.8 s. The field of view is 2".08 (EW) x 81".92 (NS), part of 
which is shown here. The pixel size is 0".16. The black region in the Stokes I map, which 
corresponds to the white region in the Stokes V map, is a pore. These images show that the 
pore lives for at least one hour. 
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Fig. 2. — Stokes I profile (top) and Stokes V profile for the region #5 (Tabled]). Wavelength 
positions from a through g define the integration ranges specified by A c , d\, d^, d% in equations 
dl]) — (j3J). a through g indicate a : X c — d 3 , b : A c — d 2 , c : A c — dx, d : A c , e : A c +di, / : A c + o? 2 , 
and g : A c + d 3 , respectively. 
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Fig. 3. — The boxes in each image indicate the region-of-interest (ROI) for the wave analysis. 
Left: SOT/SP images for the region #05 (Tabled]). The ROI with size of 1".28 (EW) x 
0".64 (NS) is located inside the pore. Right: SOT/SP images for the region #10. The ROI 
contains inter-granular magnetic Structure (IMS). The images show that the IMS is not 
maintained for 1 hour. The size of the ROI is 2". 08 (EW) x 3".52 (NS). 
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Table 1: List of observed magnetic flux concentrations. 



1 IMS: Inter-granular Magnetic Structure 

2 X-Y coordinate of the target region. X is to the West, and Y is to the North. 
3 Helio-longitudinal angle from the meridional line 
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peaks, which are all tabulated in Table [2j We analyzed 29 flux tubes, and such common 
peaks are found in 14 (48%) flux tubes, which are all tabulated in Table [TJ 

We derive the r.m.s. amplitudes of the line-of-sight fluctuation in magnetic flux 
(5$ios,rms) and velocity (Sv i os ,rms) , the line core intensity fluctuations (<5/ C ore,rms), and the 
continuum intensity fluctuations (£/ C ont,rms) at the peak periods in the power spectra. 
We also obtain phase difference between the fluctuations in the magnetic flux (4>b), the 
velocity (<fi v ), the line core intensity (0/ )CO re), and the continuum intensity (0/, C ont); 4>b — <f>v, 
4>v — 0/,core, 0/,core — 4>b, and (pi^ore — 0/,cont, all for the peak periods. The phase relations 
between the fluctuations in the magnetic flux, the velocity, and the intensity fluctuations 
are of critical importance to identify modes and properties of magneto-hydrodynamic waves 
as we will see later. 

When x n is the raw time profile (0 < n < N — 1) (N is the number of data points), 
then the complex amplitude Xk at the frequency k in the frequency domain is converted 
to the r.m.s. (root mean square) value of the wave amplitude Ak^ms and the phase 9k as 
follows: 

1 ^-4 2iiikn 

n=0 

A k , ms = y/2\X k \, (6) 

9k = arctan[ i^y ] - (7) 

We calculate these values for all the peaks, and Table [2] lists the l.o.s. magnetic flux 
^o.ios = -S 0i ios/, where -E>o,ios is the line-of-sight magnetic field and / the average filling 
factor, both of which are derived from Milne-Eddington inversion, the r.m.s. line-of-sight 
magnetic flux fluctuations (5$i os ,rms), <5< ^'° s i ' rms , the r.m.s. line-of-sight velocity fluctuations 
(^ios,rms), the r.m.s. line core and continuum intensity fluctuations normalized by the 
average intensity, corc - rms and £±S£i£Hi, an d the phase difference among magnetic, velocity, 

J core ^cont 

and intensity fluctuations; <p B ~ 4>v, <f> v ~ <t> i, core, 0/,core - 4>b, and 0/ iCor e - 4> Leant derived 
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from Eq. (j7]). 

There are 8 cases for pores and 12 cases for IMSs where magnetic, velocity, and 
intensity fluctuations have strong power at the same periods. The histograms of the phase 
difference and period for 20 such common peaks are shown in Fig. [51 The histograms for the 
phase difference show striking concentrations at around —90° for <Pb ~ <Pv an d (j) v — 4>i, C ore, 
at around 180° for <pi jCore — (fis, and at around 10° for 4>i iCor e — <fii,cont- Here, for instance, 
4>b — 4>v ~ —90° means that the velocity leads the magnetic field by a quarter of cycle. The 
periods are around 3—5 min for pores, while the periods are around 4—9 min for IMSs. 
There is no power between 134s (the detection limit due to the Nyquist criteria, see section 
2.1.) and 204s (region #04 in Table [2]). 

As pointed out in Sect. 2.2, no cross-talk should be expected in the l.o.s. magnetic 
signal from the velocity fluctuations. Furthermore, the phase difference between the 
magnetic flux and the velocity fluctuation <Pb — 4>v, if caused by the cross-talk, should be 0° 
or 180°, while the observed phase difference shows a strong concentration at around —90°. 
A similar phase relation is obtained by Bellot Rubio et al. (2000) for sunspot umbrae. On 
the other hand, Riiedi et al. (1999) and Norton et al. (1999) came to an opposite conclusion 
that the magnetic field leads the velocity by about a quarter of a cycle. 
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magnetic flux velocity intensity 




Fig. 4. — Top: Time profiles for the region #04 of Table[TJ the line-of-sight (l.o.s.) magnetic 
flux (left), the l.o.s. velocity (center), and the line core intensity (right) as defined by Eq. 
([3]). The intensity profile is normalized to the peak value of the time profile. Images of the 
region #04 are shown in Fig. [U Bottom: The power spectra of the l.o.s. magnetic flux 
(left), the l.o.s. velocity (center), and the normalized line-core intensity (right). The circles 
indicate the common, isolated peaks. 
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Fig. 5. — Left: Histograms of the phase difference between fluctuations in the magnetic 
flux, the velocity, the line-core intensity, and the continuum intensity, <$b — <p v (top, left), 
<Pv - 4>i,core (top right), 4> I<core - <p B (bottom left), and 0/ iCore - </>/, CO nt (bottom right). Solid 
lines indicate the phase difference for pores and IMSs, while dashed lines the phase difference 
for pores only. The histograms show striking concentrations at around —90° for 4>b — <Pv and 
4>v — 0/,core, at around 180° for 0/ jCO re — <Pb, and at around 10° for 0/ jCO rc — 0/,cont- Right: 
Histogram of the periods of the common peaks in the power spectra. The peak periods are 
around 3—5 minutes for pores, while the peak periods for IMSs are around 4—9 minutes. 
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Table 2: Physical parameters corresponding to the principal peak in the power spectra of all region of interests (shown 
in Table. [1]) with common peaks in the magnetic flux, the velocity, and the photometric intensity. 



r.m.s. (root mean square) l.o.s. (line of sight) magnetic flux amplitude 
2 .o.s. magnetic flux from Milne-Eddington inversion 
3 average filling factor 
4 r.m.s. l.o.s. velocity amplitude 

5 r.m.s. intensity fluctuation normalized by the average intensity for line core and continuum 
6 period 

7 phase difference between magnetic, velocity, line core intensity, and continuum intensity fluctuations 
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4. Intensity Fluctuation 



Previous authors (e.g. Bellot Rubio et al. 2000) detected fluctuations in the magnetic 
field strength and the velocity for a sunspot umbra, and obtained a phase difference of 
~ 90°. They concluded that the observed fluctuations in magnetic field strength is mainly 
caused by the opacity effect. Temperature and density fluctuations associated with the 
propagation of a hydrodynamic (acoustic) or magneto-hydrodynamic (magneto-acoustic) 
wave may cause the opacity fluctuation that moves the line formation layer upward or 
downward, resulting in an apparent magnetic field fluctuation, if the magnetic field has a 
gradient with geometrical height (dB/dz). This is called the opacity effect. 

In this section, we consider whether the observed fluctuation is due to the opacity 
effect. The photometric intensity that we observe is given by 



where T is the local temperature at the optical depth r. The intensity modulation can take 
place due either to change in the temperature or to change in the optical depth, which 
depends on the density and the temperature in the optical path. The opacity effect involves 
the second term (e~ T ). Fluctuation in intensity indicates a compressive nature of the 
fluctuation due to the first term ( y ' ) and/or to the second term (e~ T ) in eq. Thus, 
waves with low intensity fluctuation, especially those with an intensity fluctuation close to 
zero, can be considered to be a incompressible mode (such as the kink mode), while those 
with high intensity fluctuation can be considered to be a compressible mode (such as the 
sausage mode). 

The top panels of Fig. [6] show the histograms of the line core (5/ cor e,rms) and the 
continuum (5/ CO nt,rms) intensity fluctuations normalized by the average intensity 7 core and 



7 cont ; c ° rc,rms (core fluctuation) and c ° nt,rms (continuum fluctuation) for all the peaks. The 
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relation between the core and the continuum fluctuations is shown in the bottom panel 
of Fig. O The scatter plot indicates that the fluctuation at the line core is larger than 
the continuum fluctuation for all the peaks, and that the line-core and the continuum 
fluctuations are linearly correlated. A linear fitting between the line-core and the continuum 
fluctuations is given by, 

Wcont^rms = q ^W™ _ Q QQQgg (9) 
-^cont Zcore 

The cross correlation coefficient is 0.91. Fig. [5] shows that phase difference between the 
intensity fluctuation in the core and in the continuum, (0/ )CO re — 0i>ont)> has a concentration 
at around 10° ± 14°. 

We here consider the opacity effect due either to the density or to the temperature 
fluctuations. First we assume only the density fluctuation (without the temperature 
fluctuation). Magnetic field strength is smaller with height (dB/dz < 0) because of the 
canopy structure of magnetic flux tubes. Since the observations are carried out with 
25.2 to 49 degrees away from the normal, we simply assume here that the magnetic field 
strength along the line of sight decreases with height in the following discussion. The 
temperature is lower with height below the temperature minimum. When the atmosphere 
in the line formation layer is compressed (or decompressed), the line formation layer moves 
upward (downward), because the opacity along the line of sight in the flux tube increases 
(decreases). When the line formation layer moves upward (or downward), both the magnetic 
field strength and the intensity decrease (increase). Therefore, the observed magnetic field 
strength and the constant-temperature intensity fluctuation caused by the opacity effect 
should have had the phase difference of 0°, while the observed phase differences 0/, C ore — <Pb 
have a concentration at around 180°. Thus, the observed phase difference is not consistent 
with that caused by the opacity effect, if the opacity effect is caused only by the density 
fluctuation without temperature fluctuation. 
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On the other hand, the line formation layer may be compressed (or decompressed) 
under the adiabatic condition. We here consider the opacity effect due to temperature, 
assuming that the optical depth r depends only on the temperature. The dominant 
absorber in the visible wavelengths is the H~ ion (e.g. Stix, 2002). The populations of H~ 
and HI are related with the Saha equation (Rutten 1995, eq. (8.2)) 

iV(HI) , , [/(HI) 5040v , . 

log = "0.1761 - log P e + log Jjj^pj + 2-5 logT e - (10) 

where P e is the electron pressure, T e the electron temperature, x t ne ionization energy from 
H to H, iV(H~) and JV(HI) the population densities of H~ and HI, E7(H _ ) and U(EI) the 
partition function of H~ and HI. Equation ffTOl) indicates that the population of H™ depends 
highly on the temperature, and decreases with the temperature in the case of the adiabatic 
compression, while the population depends weakly on the pressure, and increases with the 
pressure in the constant temperature case. Thus, we cannot determine the population of 
H~ in the actual situation without employing a model taken into account the radiation 
exchange between the inside and the outside of the flux tubes. 

We point out that regardless of mechanism to change the opacity, the phase difference 
between the fluctuations in the magnetic field and the intensity (<fti — <$b) depends only on 
the sign of magnetic gradient along the line of sight when the line formation height moves 
upward or downward due to the lateral expansion and the contraction of the tube. The 
flux tubes that we observed were located 25.2 to 49.0 degrees away from the sun center. If 
the magnetic field strength decreases with height along the oblique line of sight, the phase 
difference between the fluctuations in the magnetic field and the intensity <pj — <pB should 
have been 0°, whereas we obtained <pi — <pB ~ 180°. Therefore, the phase relation between 
the fluctuations in the magnetic field and the intensity from the observation would not be 
consistent with that caused by the opacity effect under the assumption of the decreasing 
field strength with height along the line of sight. 
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If the effect of the adiabatic compression (or decompression; first term in equation (jSj)) 
is larger than the opacity effect due to the density and/or temperature fluctuation (second 
term in equation flS}), the phase difference between the magnetic field strength and the 
intensity fluctuation is 0° for the case of the fast-mode wave, while that is 180° for the case 
of the slow-mode wave. Thus, we can rule out the fast-mode wave, since the observed phase 
difference is close to 180°. 



In this chapter, we examine whether the observed properties of waves are consistent 
with the kink mode MHD waves (Fig. [7j). Though the magnetic and velocity fluctuations 
that we observe could be either parallel or perpendicular to the flux tubes, we here consider 
the possibility that the observed fluctuations are transverse to the magnetic field. As 
discussed in chapter 4, Fig. shows that some of the fluctuations has very small intensity 
fluctuation. Since the kink mode is essentially of non- compressive nature, those fluctuations 
with little intensity fluctuation may dominantly have properties of the kind mode. 



The dispersion relation of the kink mode neglecting gravitational stratification is given 
by (e.g. Spruit, 1981; Edwin and Roberts, 1983, Moreno-Insertis, Schussler, & Ferriz-Mas, 
1996, Ryutova & Khijakadze, 1990) 



where Ck is the phase speed of the kink mode, ui the frequency, k the wave number, va 
the Alfven speed, pi the density inside the flux tube, and p e the density outside the flux 
tube. The transverse displacement of the flux tube Sx with geometrical height z and time t 



5. 



Kink Mode MHD Waves 



5.1. 



Reflection of Kink waves 




(11) 
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Fig. 6.— Top: Histograms of continuum 5/ con t and line core 5/ C ore intensity fluctuations 
normalized by the average intensity I cont and 7 core , <57cont ' rms (left) and SIc ° rc - Lms (right) (solid 
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lines). Dashed lines indicate histograms for pores. Bottom: Scatter plot between the inten- 
sity fluctuations c . ont ' rms and c p rms . The solid line indicates the linear regression line. 
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Fig. 7. — Left:kink mode MHD wave Right: sausage mode MHD wave 
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can be expressed as Sx(z,t) = Xocos(ujt ± kz), where xq is the amplitude of the transverse 
displacement. The transverse magnetic field and velocity component are given by, 

5B = B ^^ = TB ksm(ujt±kz), (12) 
az 

d(Sx) 

Sv = — = —uo sm.(uot ± kz), (13) 
where Bq is the vertical magnetic field strength. From Eqs. ( fTTl) — ( FTBl . we obtain 

Bq iO/k 

SB = ± v / 4n(p l + p e )5v. (15) 

The phase relation of the kink mode is the same as that of the Alfven mode. Magnetic field 
is directed away from the Sun in our case. If the kink wave propagates to the direction 
same as that of magnetic field vector, minus sign should be taken, and vice versa. If a pure 
ascending or descending kink wave propagates toward the observer along magnetic field, 
phase difference between the magnetic field and the velocity fluctuations (4>b — 4>v) should, 
therefore, have been 180° or 0°, respectively. Fig. [5] shows that this is not the case. 

We then consider a superposition of the ascending kink wave and the descending 
waves, which is the reflected ascending wave at the photosphere-chromosphere boundary. 
When the ascending and the descending kink waves coexist in the line formation layer, 
the superposed wave form is determined by six variables SB U , Sv u , <p u , SB d , Svd, (fid, which 
indicate the amplitude of the magnetic field fluctuation (SB), the amplitude of the velocity 
field fluctuation (Sv), and the initial phase (cfi) of upward (subscript u) and downward 
(subscript d) waves. When magnetic field vector is toward the observer, the transverse 
magnetic field and velocity displacement of the superposed kink wave are given by 

SB = —SB U cos(ujt + (fi u ) + SB d cos(—(ut + (fid)), (16) 
Sv = Sv u cos(ut + (fi u ) + Svd cos(—(ojt + (fid))- (17) 
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Note that the phase difference between magnetic and velocity fluctuation in the ascending 
kink wave is 180°, while that in the descending kink wave is 0°. This fact is reflected in the 
sign of each term in Eqs. ffTB]) and (|17[) . We can rewrite these equations as follows: 

SB = 5B s cos{ut + (j)B), (18) 
5v = Sv s cos(ut + 4>v)i (19) 

where 8B S and 8v s are the magnetic and the velocity amplitudes of the superposed kink 
wave, and cfis and <f> v are phases of the magnetic field and the velocity of the superposed 
kink wave. In Eqs. ffTB"]) and f|T9|) . 5B S , 6v s , 4>b, and <fi v are given by: 

SB S = yfsB* + 5Bj - 2SBJB V cos(0 M - <fi d ), (20) 

8B u sm(f) u - 8B d sm<f) d 

cos0 B = — , (21) 

oB s 

. -5B u cos<j) u + 5 B d cos 4> d 

sm0B = — , (22) 

oB s 

Sv s = yj 5vl + 5vj + 25v u 5v v cos(0 u - 4> d ), (23) 

-5v u sm(f) u - 8v d sm4> d 

cos0„ = , (24) 

dv s 

. 5v u cos 4> u + 5v d cos (fi d 

sm<p v = . (25) 

dv s 

^From Eq. (fT5j) . we obtain -f^ = , 1 =. Therefore, we obtain the following relation 
among the quantities in Eqs. (1T61) and (flTj) : 

5B u ~5B d [ } 

Using Eqs. fl20|) — (1261) . the following phase difference between magnetic and velocity 
fluctuations is obtained: 

cos(0£ — (fi v ) = cos (f>B cos (fi v + sin <fi B sin <fi v 
-5BJv u + 5B d 5v d 5B u /5v u 2 2 5v u /SB u 2 2 

= Jbm, = Ibj^^ + 5Vd) = -Jb^7 { - 6B - + 6Bd) - (27) 
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This equation shows that the phase difference between the magnetic and the velocity 
fluctuations — <p v ) should be —90° or 90° when the amplitude of the reflected descending 
kink wave is exactly the same as that of ascending kink wave (i.e. 5v u = 5v d and 5B U = 5B d ). 
The observed phase relation is consistent with this prediction. 

5.2. Standing kink waves 

The transverse displacement of magnetic field line in the presence of upward (5x u ) and 
downward (6x4) kink wave is written as a function of height (z) and time (t), 

5x u (t, z) = x u0 cos(cut + kz + (j) u ), (28) 
6x d (t, z) = x do cos(ut - kz + (f) d ), (29) 

where x u o, x d o, 4> u , 4> d are the transverse amplitude and the initial phase of the magnetic field 
line fluctuation in the presence of the upward (subscript u) and the downward (subscript 
d) kink wave. When x u q = x d o = Xq, which corresponds to the case for perfect reflection, 
the transverse displacement 5x s of the magnetic field line in the presence of the superposed 
kink waves is given by 

Sx s (t, z) = 5x u (t, z) + 8x d (t, z) = 2xo cos(cut + ^ u 1" ^ d ) cos(kz + — ^ d ). (30) 

Equation (I5D1 shows that the superposed kink wave, if with perfect reflector, is a standing 
wave. Sketches of standing kink wave are shown in Fig. [HI Whether the phase difference is 
90° or —90° depends on the distance from the reflection boundary (node). 

5.3. Phase Difference 

We here give one interpretation for the concentration of the phase difference at around 
—90° (Fig. [HJ. When the ascending kink wave is reflected back at chromosphere-corona 
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boundary, and the ascending and the descending kink waves coexist in the line formation 
layer beneath the reflector, the phase difference between the magnetic and the velocity 
fluctuations should have been either 90° or —90°, while observed phase angle concentrates 
at around —90°. Whether the phase angle is 90° or —90° depends on the distance between 
the reflector and the line formation layer (Fig. [8]). The concentration at —90° indicates 
that the separation between the reflecting boundary and the line formation layer is fixed 
for all the flux tubes such that it corresponds to —90° phase difference. If we perform 
similar observations with different absorption lines with different formation height, and the 
difference in height is larger than the quarter of the wavelength (800km), this conjecture 
can be verified. 



5.4. Leakage of Poynting Flux to Corona 

The Poynting flux above the reflecting layer is the Poynting flux of the ascending kink 
wave minus the Poynting flux of the descending kink wave in the line formation layer. We 
here estimate the effective or residual upward-directed Poynting flux along a flux tube 
above the reflector. The Poynting flux of the kink wave is given by F = ^-(5B ims 8v ims ), 
so that the difference of the Poynting flux between the ascending and the descending kink 
waves is given by 

AF = — - (5B U rms^u.rms $-Bd,rms(^d,rms) , (31) 
47T 

where 5B u>Tms = 5B U /V2, 5v u , rms = 5v u /V2, 5B dtima = 5B d /V2, 6v d>imB = 5v d /y/2. Using 
Eq. ([2ZD , we can rewrite the equation as follows: 

AF = --^-{6B B tms 6v 8itms ) cos((pB - <j>v), (32) 

where 6B STTOa = 5B s /^/2, 5v S:Vms = 5v s /\^2. It turns out that the effective upward-directed 
Poynting flux is proportional to cos(0b — <p v ). 5B s ^ Tms , and Sv S}Tms in Eq. ( l32l) are related 
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to the observables, assuming that the fluctuations are transverse (i.e. normal to the flux 
tubes), 

SB s , ms = (33) 
/ sm9 

dv S!Tms = — — — , (34) 
sm V 

where 9 is helio-longitudinal angle from the meridional line. If the phase difference from 
—90° is just 6° as an exercise, i.e. <Pb — <Pv = —96°, we obtain AF = 2.7 x 10 6 erg cm -2 
s" 1 by substituting B = 1.7 x 10 3 G, 5$i os ,rms = 7.7 G, 5vi OS:Vms = 0.059 km/s, J = 0.73, 
and 9 = 29° (region #07). Therefore, even if we observe the considerable reflected wave in 
the photospheric layer with SOT/SP, there will be substantial leakage of the upward kink 
wave toward chromosphere and corona in terms of the energy flux required for the coronal 
heating (~3x 10 5 erg cm -2 s" 1 for the quiet Sun; Withbroe & Noyes, 1977). 



5.5. Seismology of Photospheric Flux Tubes 

We show in this chapter that various physical parameters that characterize the 
magnetic flux tubes are obtained simply from the amplitude and period of the magnetic 
and velocity fluctuations. We estimate the physical parameters for the region #02. The 
intensity fluctuation is 0.17—0.25% in continuum (Table [2]), and we assume that the 
observed fluctuation is due to the superposition of upward and downward kink waves. 

We define the coronal/chromospheric boundary, which is considered to be a reflector, 
to be the origin of the z-axis, which is normal to the solar surface (away from the Sun). 
A schematic behavior of the standing kink wave is shown in the left panel of Fig. [9j 
Substituting i^th, = o (without losing generality) and ^ u ~^ d = | (to make the height at 
z = the node) into Eq. ( l30l) . the transverse displacement of the flux tube is given by, 

Sx s (t, z) = 2x cos(a;t) sm(kz). (35) 
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The transverse components of the magnetic field and the velocity are given by, 

dSx 

5B s (t, z) = Bq s = 2B X()k coa(ut) cos(fcz), (36) 
05x 

Sv s (t,z) = s = —2x u>sm(uit)sm(kz), (37) 

Equations (|36|) and (1371) indicate that the phase difference between the fluctuations in the 
magnetic field and the velocity 4>b — 4>v is, 

90° for (n + §)tt < kz < (n + l)vr (sector (a) in Figl]), 

(38) 

-90° for nir < kz < (n + |)7r (sector (b) in FigEl), 
where n = —1,-2,-3,.... Equation (13"B"1) indicates that the observed phase difference 
<Pb — 4>v ~ —90° is consistent with the situation that the line-formation height is located in 
the sector (b). From Eqs. (TTTj) . (1361) . (|37|) . we have 

W w/A;. |tan(fcz)| 

/|&BJ\2|tan(^)| 2 

where |<5-B S | and |5v s | are the amplitude of the fluctuations in the magnetic field and the 
velocity, and are the function of height z. \SB S \ and \&v a \ in Eq. (flUj) are related to the 
observables, 

\bB s \ = ^ los ' rms , (41) 
/ sin 6* 

I j I V^2<5f los,rms /, n \ 

sin 6 

Assuming that the flux tubes that we observe here are in pressure equilibrium, and do not 
have a helical structure (azimuthal component), the equation for the pressure equilibrium 
for the flux tube is simply expressed as 

^ + P -k B T t = S + ^k B T e , (43) 

871 771 o7T 771 

Pe T e - Pl T l = —-(B^-B 2 e ), (44) 

07rfcn 
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where B, p, and T are the magnetic field strength, the mass density and the temperature, 
and the subscript i and e indicate the inside and the outside of the flux tube, respectively, 
m average particle mass, and ks the Boltzmann constant. From Eqs. (T4"0"j) and (PHI) , we 
can determine p; and p e assuming that outside the flux tube is field- free (B e = G), as is 
inferred by the observations. 

The line formation height in the umbra is deeper than that in the quiet Sun, because of 
the lower temperature and density (e.g. Stix, 2002). The Wilson depression for the flux tube 
with B ~ 2000 G reaches about 300 - 400 km (Deinzer, 1965; Mathew et al, 2004). The 
temperature and the average molecular weight at the height ~ —350 km is T e = 1.0 x 10 4 
K and p = 1.2 (from Table 2.4, Stix, 2002). Since the temperature inside the flux tube is 
lower than that outside the flux tube (Malt by et al, 1986), we assume = 7.0 x 10 3 K. 
We choose kz = —496° (see section 6.3 for justification to choose the value). Substituting 
m = pm p = 1.9 x 10~ 24 g, where m p is the proton mass, k^ = 1.4 x 10~ 16 erg K _1 , 
Bi = B = 1.9 x 10 3 G, 5$i os , rms = 8.8 G, 5v los , rms = 0.060 km s _1 , J = 0.75, and 6 = 29°, 
we obtain mass densities pi = 5.1 x 10~ 8 g cm -3 and p e = 2.3 x 10~ 7 g cm -3 . The 
number densities inside and outside the flux tube are n,- = — = 2.7 x 10 16 cm -3 and 

n e = ^ = 1-2 x 10 17 cm -3 , respectively. The mass density for the height of —300 400 

km is p e = 3.5 — 4.5 x 10~ 7 g cm -3 (from Table 6.1, Stix, 2002). This is consistent with our 
estimation within a factor of 2. 

We also estimate other physical parameters associated with the flux tube: (1) Alfven 
speed inside the flux tube va,% = ^fj^f i (2) plasma f3 = p '^f^ m inside the flux tube in 
the line formation layer, (3) wavelength of the kink mode L = VA,i^J~^^P, where P is 
the fluctuation period, (4) propagation time of fast magneto-acoustic wave across the flux 
tube r = , n R „ , where R is the tube radius, and (5) distance between the boundary and 
the line formation layer d = L^r. Other obvious useful parameters are the pressure scale 
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height H = , where g is the gravity in the solar surface, and the sound speed in the 

Trig 

photosphere c s = \J lk ^ T ■, where 7 is the adiabatic coefficient. Substituting R>i = 1.7 x 10 3 
G, pi = 5.1 x 1(T 8 g cm" 3 , p e = 2.3 x 1(T 7 g cm~ 3 , g = 2.7 x 10 4 cm s~ 2 , P = 312 s, 
7 = 5/3, and R = 2000 km (case # 02), we obtain v A ,i = 24 km s" 1 , (5 = 0.18, L = 3.1 x 10 3 
km, r = 75 s, d = 4.3 x 10 3 km, H = 3.9 x 10 2 km, and c s = 11 km s -1 . The propagation 
time of the fast magneto-acoustic wave across the flux tube r is less than the oscillation 
period P, and this is consistent with the assumption of the kink wave. 

Mathew et al. (2004) calculated the physical parameters (magnetic pressure, gas 
pressure, Wilson depression, and plasma f3) for a sunspot by performing an inversion to 
infrared spectro-polarimetric profiles, and derived plasma beta for the umbra (3 ~ 0.5. 
Riiedi (1992) also performed an inversion to the infrared lines, and obtained the plasma 
(3 ~ 0.25 at z = km in the plage region. The plasma beta is generally higher at z — —350 
km, following the increase in the mass density (Stix, 2002). 

As demonstrated here, we are potentially able to obtain all the physical parameters of 
the flux tube from the information on the MHD fluctuations. This indicates that seismology 
of magnetic flux tubes is possible with multiple lines corresponding to different height 
(photosphere and chromosphere) of the solar atmosphere. 

6. Sausage Mode MHD Waves 

We here consider the alternative possibility that the observed magnetic and velocity 
fluctuations are due to the longitudinal MHD waves or the slow sausage mode oscillation 
(Fig. [3 Ryutova, 2009; Defouw, 1976; Roberts and Webb, 1978; Ryutova 1981). 
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Fig. 8.— Top: A standing kink wave along magnetic field line -BO is divided into four 
parts (1) through (4) each separated by nodes and anti-nodes. Bottom: Time evolution 
of the standing kink waves. The wave evolves from (a) to (h), and goes back to (a). The 
arrows with filled box indicate velocity vector, while the arrows with circle indicate perturbed 
component of magnetic field vector. The length of the arrows indicates the magnitudes of 
the vector at certain space and time points. Schematic representation of the standing kink 
wave shows that the phase difference between magnetic and velocity fluctuations (4>b — 4>v) 
is —90° at the portions (1) and (3), and 90° at the portions (2) and (4). 
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Fig. 9. — The standing kink wave (left) and the standing slow sausage wave (right). The 
phase difference between the fluctuations in the magnetic field and the velocity (4>b ~ <Pv) is 
90° in the sector (a) and —90° in the sector (b). The arrows indicate the transverse motion 
of the magnetic fields at the anti-nodes. 
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6.1. Phase Relation for Propagating Wave 

We consider a slow mode perturbation propagating along a cylindrical flux tube, 
neglecting gravitational stratification, following Ryutova (2009). We assume that the 
magnetic and velocity fluctuations with higher intensity fluctuation (Fig. [6J may have the 
sausage-mode nature. The momentum equation perpendicular to the flux tube is given by 

B \\8B\\ 



An 



+ 5p = 0, (45) 



where the subscript means these values in unperturbed state, and 5 means perturbation 
of these values. We have the relation under the adiabatic condition 

5p = c 2 Jp, (46) 

and the flux conservation is given by 

B^8S+5B\\S = Q. (47) 
The momentum equation parallel to the flux tube is (substituting eq. (j4"5l) ) 



d5v\\ d5p 2 05 p 

p0 ^F = ~W = " c - (48) 

and the continuity equation is 

-(8pS + SSpo) + S p ^- = 0, (49) 

where S = nR 2 , Bu, p, p, vn are the cross section of the flux tube, the longitudinal magnetic 

field, the density, the pressure, and the longitudinal velocity, respectively, and c S)0 is the 
sound speed. From Eqs. ( |44l) — ( 1471) . we have 

Aixbp Atic 2 

° B W = — ^~ = ^~ ^' ( 50 ) 

-D0|| -D0|| 



SS = -S„^ = S^Sp. (51) 

B Q \\ #011 
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The continuity equation (Eq. H9]) becomes 



Taking the time derivative, and substituting Eq. ( 1481) . we have the dispersion relation, 

We therefore obtain the phase velocity of the slow sausage mode ct (c.f. Edwin and 
Roberts, 1983), 

where va is the Alfven velocity 

Hereafter we define positive as away from the solar surface. We consider a simple 
sinusoidal wave propagating upward (k > 0) or downward (k < 0) along the flux tube of 
positive (-By > 0) or negative (By < 0) polarity, 

5p = 5p cos(ut — kz) (lo = kcr), (55) 

where Sp is the amplitude of the density fluctuation. Substituting Eq. (jSTfl) . we have 

= -f-^- = -^Sp cos(ut - kz), (56) 

B \\ B \\ 



Po-^- = -c 2 s0 k5p sm(ut - kz). (57) 



and we have from Eq. (148 



Taking the integration with time (neglecting integration constant), we have 

6vn = ^^coa(wt-kz), (58* 
c T Po 

Assuming that the flux tube has an axis-symmetric sausage oscillation, the transverse 
velocity averaged over the whole pixels within the flux tube should be canceled out. Thus, 



-38- 



what we detect as a clear strong peak in the l.o.s. velocity must be longitudinal, if the 
fluctuation is due to the propagating slow sausage mode. 

^From Eqs. (|55|) . (1561) . and (158|) . we have the phase relations between the fluctuations 
in the magnetic field, the velocity, and the density, 

8p Bn , s 

— 5-, (59) 



5B\\ 47tc Sj0 



B\\ 



c 



,2 

s,0 



£=^, (61) 
dp c T p 

Equation (1601) indicates that the phase difference between the fluctuations in magnetic field 
and the velocity <$b — 4> v m the propagating wave with slow sausage mode is 0° or 180°, 
depending on the direction of magnetic field and wave propagation, whereas we observed 
4>b — 4>v ~ —90°. Thus we can rule out the possibility that the observed fluctuations are 
due to the propagating wave with slow sausage mode. 



6.2. Phase Relation for the Standing Sausage Wave 

We here consider the superposition of ascending and the descending slow sausage waves 
with the same amplitude of the density fluctuation, assuming B» > from our observation, 

5p = 5p[cos(kcTt — kz + <p u ) + cos(/ccrt + kz + <fid)] — 

op cos{kcTt H ) cos(kz H ), (b2) 

where <p u and (ftd are the initial phases of the upward and downward propagating waves with 
slow sausage mode, and k > without losing generality. From Eqs. (j5Dj) and (J6"2"|) we have 

oBu = — dp cos(kcTt H ) cos(kz H ), (63) 

Bu 2 2 
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and from Eqs. (j4"81) and fl62l) we have 

ov\\ = — —dp sm(kc T t H ) sm(kz H ). (64) 

pcx 2 2 

Equations f[6"2"j) and f[6"3"l) indicate that the phase difference between the fluctuations in the 
magnetic field and the density is 180°. 

Equations (16"3"|) and (I64p indicate that the phase difference between the fluctuations in 
the magnetic field and the velocity <Pb — <Pv is 90° or —90°, depending on the location of 
the line formation layer (right panel of Fig. [9]). The observed phase relation between the 
fluctuations in the magnetic field and the velocity is —90°. 

In the previous section, we discussed that waves with low intensity fluctuation be 
considered to be a incompressible mode (such as the kink mode), while those with high 
intensity fluctuation is considered to be a compressible mode (such as the sausage mode). 
However, Eq. ( 1621) indicates that the density fluctuation and the resultant intensity 
fluctuation are zero at the nodal points for the standing sausage wave. Thus, there may 
be cases that the standing sausage wave may not show intensity fluctuation with large 
amplitude. 



6.3. Seismology of Photospheric Flux Tubes 

We here show that the seismology of magnetic flux tubes is also possible for the sausage 
MHD oscillation. We assume that the observed fluctuation is due to the superposition of 
upward and downward compressible sausage waves for the region #05. This is justified by 
the fact that the region #05 has very high intensity fluctuation (Table [2]). 

A schematic behavior of the standing sausage wave is shown in the right panel of Fig. 
[9j Substituting i^th, = q (without losing generality) and ^ u ~^ d = | (to make the height 
at z = the node) into Eqs. fl62l — fl64l . the variations of the density, the longitudinal 
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magnetic field, and the longitudinal velocity are given by, 

5p = —5p cos(kcTt) sm(hz), (65) 
47rr 2 

SBu = —-^Sp cos(kc T t) sm(kz), (66) 
II 

c 2 

6v\\ = —^-5p smikcxt) cos(kz), (67) 

PC T 

Equations (|66|) and (1671) indicate that the phase difference between the fluctuations in the 
magnetic field and the velocity 4>b — <Pv is given by, 

90° for rm < kz < (n + h)n (sector (a) inFig. [9]), 

(68) 

-90° for (n + |)vr < kz < (n + l)vr (sector (b) inFig. EJ, 

where n = —1,-2,-3,... . Equation (1681) indicates that the observed phase difference 
<Pb ~ 4>v ~ —90° is consistent with the situation that the line-forming layer is located in the 
sector (b). 



Eqs. f l66|) and (J6"7]) are reduced to, 

47rc 2 5p| sin(/cz)|/i?|| 47rpc T | tan(kz)\ 



(69) 



l^lll c s,o^Pl cos(/cz)|/pcT B\\ 
where \5B\\ \ and are amplitudes of longitudinal fluctuations in the magnetic field and 
the velocity. B», \SBn\, and \Sv»\ are related to the observables, 

B\\ = B , (70) 

/ cos# 

cos 9 



Since c s = A / 7fcflT and va ' 



c y A ik B TB\ 



c 2 s + v\ Airp^k B T + B\m 



(73) 
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^From Eqs. (1691) and (1751) . we have 



/ Jtffijll x ^ = (4vrp) 2 7 fc i? r|tan(^)| 2 ^ 
V|<Ju[||/ 47Tfnk B T + B\ m ' 1 ] 



Equation (1741) leads to a second order equation for p, 



dip 2 - a 2 p - a 3 — 0, (75) 

ax = (47r|^|||) 2 o'fc i3 T| tan(^)| 2 , (76) 

a 2 = 47T7fc B T|5fi|| | 2 (77) 

a 3 = (Si||<JS|||) 2 m (78) 



Since p > 0, we can take only p = a2+ V^2+ 4 _^±_ This indicates that we can determine 
the mass density inside the flux tube with the additional knowledge of tan(kz) for the 
line-forming height. However, there are multiple solutions due to ambiguity in tan(fcz). 
The region that we chose for the photospheric seismology (region #02 with the assumption 
of the kink wave and #05 with the assumption of the sausage wave) are both pores, whose 
magnetic field strength is almost the same. We assume that the parameters of the flux tube 
(the mass density, plasma beta, and Alfven velocity) and distance between the boundary 
and the line formation layer derived from the analysis of the region #05 (sausage-wave 
dominant) should be consistent with those derived from the analysis of the region #02 
(kink-wave dominant, section 5.5). The choice of kz = —619° for Eq. (T40l) and kz = —496° 
for Eq. (ITS"]) in the following exercise is based on the assumption. 



Substituting B = 1.7 x 10 3 G, <J$ ]oB , riM = 13.8 G, <k>i os>rms = 0.12 km s" 1 , / = 0.81, 
6 = 29°, 7 = 5/3, k B = 1.4 x 10~ 16 erg K _1 , T = 1.0 X 10 4 K, and m = 1.9 x 10~ 24 g, 
we obtain the mass density inside the flux tube p = 0.8 x 10~ 7 g cm -3 . We then derive 
the values associated with the flux tube; (1) Alfven speed va = (2) plasma beta 

P = P B ^/%™ ; (3) phase speed of the slow sausage mode Ct, (4) wavelength of the slow 
sausage mode L = ctP, where P is the observed oscillation period, and (5) distance between 
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the boundary and the line formation layer d = Substituting P = 294 s, we obtain 

va = 23 km s -1 , j3 = 0.18, ct = 8.5 km s™ 1 (sound speed c s = 11 km s™ 1 ), L = 2.5 x 10 3 
km, and d = 4.3 X 10 3 km. 

The set of parameters derived here satisfy the condition that p (or p^), (3, va, and d 
derived here are consistent with those derived in section 5.5. The distance between the 
boundary and the line formation layer d = 3.6 x 10 3 km is consistent with the distance 
between the line-formation height and the transition region. This indeed indicates that the 
transition region is the reflecting layer for such waves. Note that the wavelength L is much 
larger than the scale height H = ~ 3.9 x 10 2 km, and the effect of the gravity has to 
be taken into account for more rigorous treatment. 



7. Discussions 

We have detected clear signatures of the MHD waves propagating along the magnetic 
flux tubes in a form of velocity, magnetic and intensity sinusoidal waves with exactly the 
same period. One or two strong and sharp peaks with common periods in the power spectra 
of the l.o.s. magnetic flux, the l.o.s. velocity, and the intensity time profiles are evident in 
the pores (8 peaks) and the IMSs (12 peaks). We note that only about half of the observed 
flux tubes have such common peaks. Periods of the peaks concentrate at around 3—6 
minutes for pores and 4—9 minutes for IMSs. Phase difference between the l.o.s. magnetic 
flux (4>b), the l.o.s. velocity (<f> v ), the line core intensity (0/ jC ore), and the continuum 
intensity (<f>i j00n t) have striking concentrations at around —90° for <p B — <p v and (fi v — 4>i, COTe , 
around 180° for 0/ ;CO rc — 4>b, and around 10° for 4>i yCOVe — 0/ iC ont (Fig. EJ). These fluctuations 
are associated with the intensity fluctuations <5/c ° nt - rms and <5Jc ° rc - rms . The amplitude of the 

Jcont i core 

intensity fluctuations amount to 0.1— 1.0 % of the average intensity level. Some flux tubes 
have a very small intensity fluctuation, and the wave mode for such flux tubes is considered 
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to be the incompressible kink mode. On the other hand, flux tubes with higher intensity 
fluctuation may have the compressible sausage mode. 

The phase relation 4>i — 4>b ~ 180° from the observation would not be consistent 
with that caused by the opacity effect (e.g. Bellot Rubio et al., 2000), if the magnetic 
field strength decreases along the line of sight toward the observer. We propose that 
the longitudinal and/or transverse MHD waves propagating along the flux tube are 
responsible for the fluctuations. The observed phase difference 4>b — 4>v ~ —90° is 
consistent with the phase relation of the superposition of the ascending and the descending 
kink wave. This indicates that the ascending kink wave is substantially reflected at the 
chromospheric-coronal boundary. The superposed waves have the property of the standing 
waves. In addition to the standing kink mode, the observed phase relation between the 
fluctuations in the magnetic flux and the velocity is consistent with the phase relation for 
the superposition of the ascending and the reflected descending slow-mode sausage waves. 

So far our analysis is based on the assumption that the either the kink mode or the 
sausage mode is dominant in the flux tubes. Both the kink mode and the slow sausage 
mode may be excited in the same flux tube. Torsional waves are not discussed in this paper. 
The region of interest encompasses the entire magnetic flux concentrations in the spatial 
and temporal domain (i.e. in the case of IMSs), and we average the physical parameters 
inside the ROI. Thus, we are probably unable to detect the torsional Alfven waves, even if 
they exist, because the perturbation of the magnetic flux and the velocity is averaged over 
the whole flux tube, and are canceled out. 

We derive the various physical parameters of the flux tubes only from the observed 
period and the amplitudes of magnetic and velocity oscillations. Such parameters include 
(1) mass density inside and outside the flux tube, (2) plasma (3 inside the flux tube, (3) 
Alfven speed inside the flux tube, (4) phase speed, (5) wavelength, (6) distance between the 
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boundary and the line formation layer, and (7) propagation time of fast magneto-acoustic 
wave across the flux tube. In the examples presented in this paper, we choose similar 
sets of kz as defined in sections 5.5 and 6.3 for both cases (the kink-wave dominant case 
and the sausage-wave dominant case) such that the derived physical parameters of the 
flux tubes coincide. The choice determines the distance d between the boundary (node) 
and the line formation layer. The derived mass density outside the flux tube is consistent 
with that of the standard solar model in the case of the kink wave. Note that we can 
not derive the mass density outside the flux tube in the case of the slow sausage mode, 
because the flux tube is not in the pressure equilibrium. This exercise demonstrates that 
the seismology of magnetic flux tubes is possible with the observations of the oscillation 
period and amplitudes for various photospheric and chromosheric lines, and may open a 
new channel for the diagnostics of the magnetic flux tubes. 

Magnetohydrodymanic waves are believed to play a vital role in the acceleration and 
heating of the fast solar wind. However, it has been thought that the Alfven speed rapidly 
increases with height due to the rapid decrease in the plasma density, resulting in significant 
reflection at the chromosphere- corona boundary. We indeed show that this may be the case 
in this paper: the upward propagating kink and/or sausage waves must be significantly 
reflected back above the line formation layer. Deviation in the phase difference between 
the magnetic and velocity fluctuations from —90° as seen in Fig. 6 may indicate residual 
waves propagating to the corona. Indeed, the upward Poynting flux above the reflecting 
layer is estimated to be 2.7 x 10 6 erg cm" 2 s -1 in one case (kink wave), and is by no means 
negligible flux in terms of heating and acceleration of the upper atmosphere. 

Tsuneta et al (2008b) conjectures that a rapid decrease in the magnetic field strength 
associated with the rapidly expanding flux tube near the chromosphere-corona boundary for 
the polar kG patches reduces the vertical change in Alfven speed, and the Alfvenic cutoff 
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frequency be lower in the polar flux tubes. Magnetohydrodymanic waves generated in the 
photosphere may be more efficiently propagated to the corona through such fanning-out 
flux tubes with large expansion factor observed in the polar coronal holes. On the other 
hand, the observations presented here suggest significant reflected waves. It would therefore 
be interesting to see whether the reflectivity of the magnetohydrodymanic waves depends 
on the locations or environment e.g. coronal holes vs the quiet Sun. 

Two interpretations addressed here (kink and sausage MHD modes) cannot be 
distinguished in the present study. The flux tubes that we analyzed are located with angular 
distance of 25—49° from the Sun center for high sensitivity magnetic observations. It is 
important to compare the wave properties for the flux tubes located further away from the 
Sun center with those of the flux tubes around the disk center to separate individual modes 
of waves (Norton et al. 2001). These topics will be addressed in our subsequent paper. 
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Appendix 

We here determine the conversion coefficient A to convert the circular polarization CP 
to the l.o.s. magnetic flux in equation (j3j). Fig. [TU1 shows the scatter plot for the circular 
polarization derived by equation (pQ) with the l.o.s. magnetic flux. The l.o.s. magnetic flux 
$o,ios = -B 0) ios/ is determined from the line-of-sight magnetic field strength B j os and the 
filling factor /, both obtained with the Milne-Eddington inversion. The data used here is 
the plage region ^05 (Tabled]) taken at 19:15 UT on 2007 Feb. 3. We notice a good linear 
correlation between the two quantities. The linear regression is given by 

CP = (4.16 x 10~ 5 )$ ,ios + 0.0016. (79) 

The correlation coefficient is 0.96. We use the conversion coefficient A = 4.16 x 10~ 5 G _1 
for the analysis presented in this paper. 
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Fig. 10. — Scatter plot indicating the relation between the l.o.s. magnetic flux $o,ios ( see 
text) and the circular polarization CP as defined in equation (CQ). Solid line indicates a linear 
regression line. 
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